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We investigate the effect of a coupling between dark matter and dark energy on the rates for 
the direct detection of dark matter. The magnitude of the effect depends on the strength k of this 
new interaction relative to gravity. The resulting isothermal velocity distribution for dark matter 
in galaxy halos is still Maxwell-Boltzmann (M-B), but the characteristic velocity and the escape 
velocity are increased by y/1 + k'^. We adopt a phenomenological approach and consider values 
- - ■ of K near unity. For such values we find that: (i) The (time averaged) event rate increases for 

Jir ' light WIMPs, while it is somewhat reduced for WIMP masses larger than 100 GeV. (ii) The time 

I dependence of the rate arising from the modulation amplitude is decreased compared to the standard 

. M-B velocity distribution, (iii) The average and maximum WIMP energy increase proportionally 

' to 1 + K'^, which, for sufficiently massive WIMPs, allows the possibility of designing experiments 

^1^' measuring 7 rays following nuclear de-excitation. 
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The combined MAXIMA-1 |1|, BOOMERANG 0, DASI i3|] and COBE/DMR Cosmic Microwave 
Background (CMB) observations imply that the Universe is flat |^ and most of its energy content 
^\ is exotic 0). These results have been confirmed and improved by the recent WMAP data 0- The 
1 deduced cosmological expansion is consistent with the luminosity distance as a function of redshift 
of distant supernovae According to the scenario favored by the observations there are 

various contributions to the energy content of our Universe. The most accessible energy component 
is baryonic matter, which accounts for ~ 5% of the total energy density. A component that has not 
been directly observed is cold dark matter (CDM)): a pressureless fluid that is responsible for the 
O ^' growth of cosmological perturbations through gravitational instability. Its contribution to the total 
O . energy density is estimated at ^ 25%. The dark matter is expected to become more abundant in 
' extensive halos, that stretch up to 100-200 kpc from the center of galaxies. The component with the 
biggest contribution to the energy density has an equation of state similar to that of a cosmological 
constant and is characterized as dark energy. The ratio w — p/ p is negative and close to —1. This 
, component is responsible for ~ 70% of the total energy density and induces the observed acceleration 
■ of the Universe The total energy density of our Universe is believed to take the critical 

value consistent with spatial flatness. 

Since a non-exotic component cannot exceed 40% of the CDM 11], there is room for a compo- 
nent consisting of exotic weakly interacting massive particles (WIMPs). Supersymmetry naturally 
provides candidates for these dark matter constituents 0, 0| . In the most favored scenario of su- 
persymmetry, the lightest supersymmetric particle (LSP) can be described as a Mai or ana fermion, 
a linear combination of the neutral components of the gauginos and higgsinos |l2l|-|l4|. In most 
calculations the neutralino is assumed to be primarily a gaugino, usually a bino. Even though there 
exists firm indirect evidence for a halo of dark matter in galaxies from the observed rotational curves, 
it is essential to detect such matter directly Until dark matter is actually detected, we 

will not be able to exclude the possibility that the rotation curves result from a modification of the 
laws of nature as we currently view them. The direct detection will also reveal the nature of the 
constituents of dark matter. 

The possibility of direct detection, however, depends on the nature of the dark matter constituents. 
Since the WIMPs are expected to be very massive {mwi mp ^ 30 GeV) and extremely non relativistic 
with average kinetic energy (T) ~ 50 KeV (mvi//Afp/100 GeV), they are not likely to excite the 



nucleus. As a result, they can be directly detected mainly via the recoiling of a nucleus {A, Z) in 
elastic scattering. The event rate for such a process can be computed from the following ingredients: 



1. An effective Lagrangian at the elementary particle (quark) level obtained in the framework of 
the prevailing particle theory. For supersymmetry this is achieved as described in refs. iHili, 
for example. 

2. A well defined procedure for transforming the amplitude obtained using the previous effective 
Lagrangian from the quark to the nucleon level, i.e. a quark model for the nucleon. This step 
in SUSY models is non-trivial, since the obtained results depend crucially on the content of 
the nucleon in quarks other than u and d. 

3. Knowledge of the relevant nuclear matrix elements |l7lll8j . obtained with reliable many-body 
nuclear wave functions. Fortunately, in the case of the scalar coupling, which is viewed as the 
most important, the situation is a bit simpler, as only the nuclear form factor is needed. 

4. Knowledge of the WIMP density in our vicinity and its velocity distribution. Since the essential 
input here comes from the rotational curves, dark matter candidates other than the LSP 
(neutralino) are also characterized by similar parameters. 

In the past various velocity distributions have been considered for the dark matter gas in our 
galaxy. The most popular one is the isothermal Maxwell-Boltzmann (M-B) velocity distribution 
with (v^) = 3v^ ~ 3^o/2i where v"^ = (w^) — {Vy) = {vD and vq is the velocity of the sun around 
the galaxy, i.e. uq — 220 km/s. Extensions of the M-B distribution have also been considered, in 
particular these that are axially symmetric with enhanced dispersion in the galactocentric direction 
|l9t l2n||. In such distributions an upper cutoff v^sc — 2.84 wq is introduced by hand, in order to 
eliminate velocities above the escape velocity. 

Non-isothermal models have also been considered. Am ong these one should mention the ones with 
late infall of dark matter into the galaxy, i.e caustic rings ^2S|-'23|, dark matter orbiting the Sun 26], 
Sagittarius dark matter j^^] . The velocity distribution has also been obtained in " adiabatic" models 
employing the Eddington proposal l28l|-|8lj. In such an approach, given the density of matter, one 
can obtain a mass distribution that depends both on the velocity and the gravitational potential. 
Evaluating this distribution in a given point in space, e.g. in our vicinity, one obtains the velocity 
distribution at that point in a self-consistent manner. Unfortunately this approach is applicable 
only if the density of matter is spherically symmetric. 

In the present work we will consider another variant of the isothermal M-B distribution, that 
results when the dark matter interacts with the dark energy I s^. l33l| . The difficulty with explaining 
the very small value of the cosmological constant that could induce the present acceleration has 
motivated the suggestion that this energy component is time dependent [sS . lisf . In the simplest 
realization, it is connected to a scalar field </> with a very flat potential. The vacuum energy associated 
with this field is the dark energy that drives the acceleration. If such a field affects the cosmological 
evolution today, its effective mass must be of the order of the Hubble scale, or smaller. 

It is conceivable that there is a coupling between dark matter and the field responsible for the 
dark energy |36l |. In such a scenario it may be possible to resolve the coincidence problem, i.e. the 
reason behind the comparable present contributions from the dark matter and the dark energy to 
the total energy density. The presence of an interaction between dark matter and the scalar field 
responsible for the dark energy has consequences that are potentially observable. The cosmological 
implications depend on the form of the coupling, as well as on the potential of the field [S^. If 
the scale for the field mass is set by the present value of the Hubble parameter, then the field is 
effectively massless at distances of the order of the galactic scale. Its coupling to the dark matter 
particles results in a long range force that can affect the details of structure formation |8^-[4ll|. 

The attraction between dark matter particles mediated by the scalar field is expected to modify 
the distribution and velocity of dark matter particles in halos, with implications for dark matter 
searches. A careful analysis indicates that the distribution remains Maxwell-Boltzmann, but with a 
potentially larger characteristic velocity [s^lssf . This has two consequences: 



• The total detection rate is reduced for large WIMP masses (above 100 GeV). This occurs 
because the velocity distribution is shifted to higher values. As a result, such a distribution 
tends to favor a high energy transfer to the nucleus. The nuclear form factor tends to suppress 
the high energy transfer components, resulting in an overall suppression. 

• The modulation effect, i.e. the periodic dependence of the rate on the Earth's motion, is 
reduced. This is unfortunate, because the modulation is viewed as a good signature against 
the background. 

• As the average WIMP velocity increases, the average WIMP energy increases as well. The 
kinetic energy becomes 

m.50a + .^,KeV^. ,1, 

Thus, for m\YiMP = 200 GeV, one finds (T) ~ 0.32 and 1.3 MeV for = 1 and 3 respectively. 
Since a value ~ 1 cannot be excluded from the available constraints, there is, in this 
case, a reasonable possibility for exciting the nucleus. In such a scenario the previous two 
disadvantages are not relevant, as they are connected with nuclear recoil experiments. This 
possibility is indeed good news, because measuring the de-excitation 7 rays is a much simpler 
task than the detection of the recoiling nuclei. 

The above conclusions depend only on the velocity distribution and nuclear structure and are inde- 
pendent of the specific nature of the WIMP. 

It must be emphasized that it is not easy to construct extensions of the Standard Model that 
include a dark energy field coupled to dark matter. The main obstruction is related to the necessity 
to keep the mass of the field of the order of the present Hubble scale after radiative corrections. A 
large coupling to the dark matter field induces significant loop corrections to the potential of the 
dark energy field, resulting in a large mass the other hand, it is reasonable to expect that 

the resolution of the coincidence problem will require a coupling that is not much smaller than the 
gravitational one. For this reason our analysis will be essentially phenomenological. We will assume 
that the dark energy field has a mass of the order of the Hubble scale and a coupling to the dark 
matter of gravitational strength. Explicit models that realize these assumptions are given in refs. 

Eilii. 



INTERACTION BETWEEN DARK MATTER AND DARK ENERGY 

We consider an interaction between the scalar field and the dark matter particles that can be 
modeled through a field-dependent particle mass. The action takes the form 



d^x^ ( M^i?- i^'^-a^^^,^- C/(0)^ m{c^{x,))dn, (2) 



2' 



where dri = \/ —g^u{xi)dx'ldx^ and the second integral is taken over particle trajectories. Variation 
of the action with respect to results in the equation of motion 

[V^ d.^) - — T';,, (3) 



where the energy-momentum tensor associated with the gas of particles is 

T^'-^-^V / dn m(0(x.)) ^^6^^\x-x,). (4) 

We are interested in static spherically symmetric configurations, with the scalar field varying 
slowly with the radial distance r. Our treatment is relevant up to a distance ri ~ 100 kpc beyond 



which the dark matter becomes very dilute. For r ^ ri we expect that (j> quickly becomes constant 
with a value close to (j>{ri) = 4>i. This is the value that drives the present cosmological expansion. 
Here we assume that the cosmological evolution of 0i is negligible for the time scales of interest, so 
that the asymptotic configuration is static to a good approximation. 

We approximate: to((/)) ~ r7i((/)o) + [d'm{(j)o) / d(p] 54> = mo + TOq S4>, with 0o the value of the field 
at the center of the galaxy (r = 0). We work within the leading order in Scf) and assume that 
m' /m ~ mQ/mo for all r. Also dU/dcj} can be approximated by a constant between r = and 
r ~ oo. For the scalar field to provide a resolution of the coincidence problem, the two terms in 
the r.h.s. of Eq. JSJ must be of similar magnitude in the cosmological solution. This means that 
dU/d(j) must be comparable to (m.o/TOo)poo- We expect poo to be a fraction of the critical density, 
i.e. Poo ~ 3 keV/cm'^. On the other hand, the energy density in the central region of the static 
solution (r ^ 100 kpc) is that of the galaxy halo 0.4 GeV/cm'^ for our neighborhood of the Milky 
Way). This makes dU/d4) negligible in the r.h.s. of Eq. ^ for a static configuration. The potential 
is expected to become important only for r — > cx), where the static solution must be replaced by the 
cosmological one. Similar arguments indicate that we can neglect U relative to p. Also the scalar 
field must be effectively massless at the galactic scale. For these reasons we expect that the form of 
the potential plays a negligible role at the galactic level. Our analysis can be carried out with U — Q 
and is model independent. 

We treat the dark matter as a weakly interacting, dilute gas. We are motivated by the phe- 
nomenological success of the isothermal sphere j4S] in describing the flat part of the rotation curves. 
We do not address the question of the density profile in the inner part of the galaxies (r ^ 5 kpc). 
We approximate the energy-momentum tensor of the dark matter as T'^ = diag(— p,p,p,p) with 
p{r) = p{r) (tj^) = m((j>{r)) n{r) {v^). The dispersion of the dark matter velocity is assumed to be 
constant and small: (u^) ^ 1. The gravitational field is considered in the Newtonian approximation: 
goo — 1 + 2<i>, with $ = O (mQ(5(/)/mo). In the weak field limit and for p p, the conservation of 
the energy-momentum tensor gives 

P' = -P1>'-P^W, (5) 
mo 

with the prime on p, 6(j) denoting a derivative with respect to r. Integration of this equation gives 
n ~ noexp - im'o/mo)S(j)/ (vj)) . 

With the above assumptions we obtain the equations of motion 

2 1 

+ r " 4M2 P° ^""P ~ ^^'^^ ' 

and 

2 m' 

{S(t>)" + - {5(t>y = ^po exp - aS^) , (7) 

r mo 

where M — (IQ-kGn)^^^^ is the reduced Planck mass, po = mono the energy density of dark matter 
at r = 0, a l/{v'^), and a — m'o/ {mo{v'^))- We emphasize that, even though |$| <C 1, the 
combination *I'/(f^), that appears in the exponent in the expression for the number density n, can 
be large. Similarly, the expansion of the mass around the value mo = m((j)o) assumes the smallness 
of the dimensionless parameter \m'oS(j)/mo\. However, the combination ad4> = (mg(50/mo)/(t'^), that 
appears in the exponent, can be large. 



A linear combination of Eqs. gives 



d'^u 2 du 

+ -— -Hexpw, (8) 
dz^ z dz 



where u — — — a(50, z = (3r and f3'^ ~ {1 + K^)apo/4M^. The parameter 



k2 = 47Vf2 („j^/„jg)2 



(9) 



determines the strength of the new interaction relative to gravity. The solutions that are regular for 
small z approach the form 



In 
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d\ cos -— In z + ^2 sin 




(10) 



for large z. Another linear combination of Eqs. 10, Q gives 

dz^ z dz ' 

with w = — K^a$ + aS(j). The solution of this equation is w = cq + ci/z. 

The velocity w of a massive baryonic object in orbit around the galaxy, at a distance r from its 
center, can be expressed as 



V \^ r^' z f du dw 



2 \ dz dz 



(12) 



where 



..2 ^ /„.2 



2 



1 + 



(vi). (13) 



The asymptotic form of u{z), w{z) indicates that w ~ for large z. The dominant correction to the 
leading behavior arises from the term ~ l/-\/z in Eq. (|10|l . The function v{z) gives a higher order 
correction. This simple analysis indicates that the approximately flat rotation curves outside the 
galaxy cores are a persistent feature even if the dark matter is coupled to a scalar field through its 
mass. If the new interaction is universal for ordinary and dark matter, the experimental constraints 
impose <C 1. In this case, it is reasonable to expect a negligible effect in the distribution of matter 
in galaxy halos. However, if interacts only with dark matter, as we assume here, this bound can 
be relaxed significantly. 

A massive particle in orbit around the galaxy, at a large distance r from its center, has a velocity 
given by Eq. (|13|l . We can use this expression in order to fix (w^) for a given value of k. The effect of 
the new scalar interaction is encoded in the factor k^. When this is small, the velocity of an object 
orbiting the galaxy is of the order of the square root of the dispersion of the dark matter velocity. 
If is large, the rotation velocity can become much smaller than the typical dark matter velocity. 

The allowed range of k is limited by the observable implications of the model that describes the 
dark sector. It is reasonable to expect that the resolution of the coincidence problem through an 
interaction between dark matter and dark energy will have to rely on a coupling not significantly 
weaker that gravity. It seems unlikely that a coupling <C 1 can lead to a cosmological evolution 
drastically different from that in the decoupled case. 

The dependence of the mass of dark matter particles on an evolving scalar field during the cosmo- 
logical evolution since the decoupling is reflected in the microwave background. The magnitude of 
the effect is strongly model dependent. In the models of ref. 37, 43] the observations result in the 
constraint ^ 0.01. In the model of ref. |4lj the scalar interaction among dark matter particles 
is screened by an additional relativistic dark matter species. As a result, the model is viable even 
for couplings ~ 1. A similar mechanism is employed in ref. this model the interaction 

between dark matter and dark energy becomes important only during the recent evolution of the 
Universe. In general, an interaction that is effective for redshifts z ^ 1 — 2 is not strongly constrained 
by the observations. 

Independently of the value of k^, the interaction of dark matter with the scalar field associated 
with dark energy does not destroy the approximately fiat profile of the rotation curves. Other 
considerations, however, could constrain the coupling k^. The dispersion of the dark matter velocity 
is (v^) = (1 + K^)w^/2. For a value of Vc deduced from observations, {v^) increases with k. For 
sufficiently large k, it seems possible that Vc may exceed the escape velocity from the galaxy. It 
turns out, however, that this is not the case. Outside the core of the galaxy and for r ^ ri, the 



binding potential for a dark matter particle is $ + {a/a)6(j). For large r, Eq. implies that 

V = — K^a$ + aS(j) = constant. The binding potential becomes (1 + k^)<1> = (1 + ln(r/ri), 
where we have omitted an overall constant. For a particle at a distance from the center of the 
galaxy, the escape velocity becomes 

vl, = 2{l + K')[\n{n/r,) + l]. (14) 

The value of Vesc is larger than the standard one by a factor (1 + k^), so that (v^) remains 
substantially smaller than Wg^^ for r, <C ri. A particle that does not interact with the scalar field is 
bound only by the the potential <&. However, the scale of its velocity is set by Vc, so that again it 
cannot escape. 



THE VELOCITY DISTRIBUTION OF DARK MATTER 



In the previous section we saw that in isothermal models the dark matter velocity distribution 
with respect to the galactic center is M-B: 

/(«) = ^^cxpf-4), (15) 

where = = (1 + k'^)v1, with Vc the observed rotation velocity of a baryonic object in orbit 
around the galaxy. This means that the dispersion of the dark matter velocity is proportional to 
1 + K^, where k is the coupling between dark matter and dark energy, given by Eq. We have 
also assumed that the ordinary baryonic matter does not couple to the dark energy field and is not 
affected by its presence. We impose an upper bound Vf, on the dark matter velocity, equal to the 
escape velocity given by Eq. (fT^ . We express it as 

Vb^nVescfl, (16) 

where Vesc,o is the escape velocity for k = and = 1 + k^. 
In the local frame this velocity distribution takes the form 

/(^) = ^^— cxp 2 , (17) 

where 

ve = i}q + t;i(sinQ! x — cosa C0S7 y + cosa sin7 z ). (18) 

We have chosen a coordinate system in which the polar z-axis is along the direction of motion of the 
Sun, the x-axis is radially out of the galaxy and y = z x x. The velocity of the Earth is Ve- The 
velocity of the Sun around the center of the galaxy is vq z (with v^f ~ 220 km/s). The magnitude of 
the velocity of the Earth relative to the Sun is vi ~ 30 km/s. The quantity 7 ~ 7r/6 describes the 
orientation of the ecliptic with respect to the galactic plane. (The angle between the normals to the 
two planes is 7r/2 — 7 ~ 7r/3.) The parameter a denotes the phase of the Earth, (a = on June 
2nd.) 

In the standard scenario, in which there is no interaction between dark matter and dark energy, we 
have Vm = and Vh = ijescVa with y^sc — 2.84. In the scenario we are considering both parameters 
are scaled up by the same factor, i.e. Vm — nvg and Vb = n y^sc vo, with n — + The standard 
M-B distribution is a special case of our model with n = 1. In the present work we treat n as a free 
parameter, which we do not expect to be much larger than unity. We will consider values as large 
as n = 2 (that corresponds to k = -y/S) and study the implications for direct dark matter detection. 
The distribution function can be written as: 

I I -y2-2((25)2+a\/l-52 coa(0) =10(0)25 + 003(0) ( (y£ + 2) 3i„(T) - bv'T^P' <!°a(T) sin( *) ) 2i + Hj + 1) 

f{y,i,(l),S,n) = — i=^e ^ , 

(19) 



where is the azimuthal angle, ^ the cosine of the angle between v and Vq, y = v/vq, S 
sin7Ui/uo ~ (l/2)(30/220) = 0.068. The integral over can be done analytically to yield: 

2 I -B^-2((2a)2+coaQ((H5 + 2)ainT)2« + y5 + l) 



/(y,e,<5,n)- -^^e ^ ^ /o [2Sy^il-e)il-cos^asm^-f)j , (20) 

where Io{x) is the well known modified Bessel function. The various variables are constrained by: 



+ 2 ((2(5)2 + cos a{y^ + 2)S + y^ + I) < ny^sc- (2f ) 

From the kinematics of the WIMP-nucleus collision we find that the momentum transfer to the 
nucleus is given by 

q = 2/irW cos6', (22) 

where 9 is the angle between the WIMP velocity and the momentum of the outgoing nucleus, and fir 
the reduced mass of the system. Instead of the angle 6 one can introduce the energy Q transferred 
to the nucleus, Q = / {2Amp) (Anip is the nuclear mass). Thus 



A 



Til 

2 sin 6* cos 6'd6' = 'f-^dQ. 

Furthermore, for a given energy transfer the velocity v is constrained to be 



/ QAnip 1 

V > Vmin , Vrmn = \ ' • (23) 

\ 2 fir 

We will find it convenient to introduce, instead of the energy transfer, the dimensionless quantity u 

u . i(,5)^ - I; ' = ^ = 4.1 X 10^ keV, (24) 

where b is the nuclear (harmonic oscillator) size parameter. 

It is clear that for a given energy transfer the velocity is restricted from below. We have already 
mentioned that the velocity is bounded from above by the escape velocity. We thus get 



aVu < y < nye 



V2firbvQ\ , (25) 



2sme cos 9d9 = -dy. (26) 

y 



THE DIRECT DETECTION EVENT RATE 

The event rate for the coherent WIMP-nucleus elastic scattering is given by [23I li^ IItI li^ : 



p(0) m 



R = — — ^fcoh{A,fir{A))al^o (27) 



with 



^ ,,,, lOOGeV 

fcoh{A, fir{A)) = 



^ir{A) 
flrip) 



1 



A tcoh (1 + hcohcosa) (28) 



In the above expression o-p^o is the WIMP-nucleon scalar cross section, p(0) the WIMP density in 
our vicinity, m^o the WIMP mass , m the target mass, A the number of nucleons in the nucleus 



and {v'^) = 3wq/2 the average value of the square of the WIMP velocity for n = 1. The number of 
events in time t is: 



The quantity of interest to us is r = tcoh (1 + hcoh cos a), which contains all the information regarding 
the WIMP velocity distribution and the structure of the nucleus. It also depends on the reduced 
mass of the system. 

The event rate is proportional to the WIMP flux, i.e. proportional to the WIMP velocity. In 
Eq. (|29|l we have chosen to normalize the event rate using the velocity dispersion for n — I, i.e. 
a/ (f 2) = ^,/3i/2yp. As a result a compensating factor of v/ ^ (w^) is included in r. It is not difficult 
to show 0, El EHH that 

pyesc 2 rio{y,a) 

P'i^) XoV^y'dy fiy,^,6,n)d^, (30) 



a. 



s 



du Ja^ V 3 2/2 

where F{u) is the nuclear form factor. In the integrand we have displayed explicitly all the factors 
of y in order to keep track of their origin. The first one comes from the flux, the second from 
the transformation (|26|l and the last is the usual phase-space factor. The quantity Co(2/jQ^) enters 
because in some region of the velocity space the upper value of ^ is restricted so that the condition 
(|21|l is satisfied. The above expression can be cast in the form: 



(31) 



We have seen that the parameter a depends on the nucleus, wq and the WIMP mass. 

In spite of the complications arising from the condition 1)21(1 . by taking the leading order of the 
modified Bessel function in Eq. H2U|I . we were able to get an analytic expression for "^{x^a) as 
follows: 

*(a;,a) = (1 - 6(2: - Uesc + 1 + (^cosa))\E'<(x, a) + 6(2; - yesc + 1 + 5 cosa)^y{x,a) (32) 

2[d cos a + 1) 

\ Prf I ~yea<: + ^ COS O+l 

f 2{dcosa + l) -yesc\ \ " 



erf 



\ n J 2((5cosa + l) 



2e 



Vesc 

(a:-y„c) _ f -y„c+i? cos g+l \ , ^^j- cos g + l ^ 



2(0 cos a + 1) 

where x is a short hand notation for a^/u and 8(x) is the Heavyside step function. 

By performing a Fourier analysis of the function "^{x,a), which is a periodic function of a, and 
keeping the dominant terms we find: 



dr [2 

— = \ —a^F'^{u) \ {a\/u) + H{a^/u) cos a + H2(a^/u) cos 20;] 
du \ 3 



(35) 



Sometimes we will consider separately each term in the above expression by writing: 
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Figure 1: The function ^0(2;) as defined in the text. From left to right n = 1 and 2. The area under the 
curve is roughly independent of n. Due to the nuclear form factor, not all the range of u is exploitable in 
direct WIMP detection. For ^^^/there is effectively a cut off value indicated by a dotted line, a fine line and 
a thick line for a WIMP mass of 30, 100 and 200 GeV respectively. The exploitable area under the curve 
decreases as n increases. 



Before proceeding further by considering a special target, it is instructive to concentrate on the 
dependence of ^0(2;) and the modulation H{x) on the parameter n of the M-B distribution. For 
this purpose we exhibit the function 'I'o(a;) in Fig. ^ the function H{x) in Fig. B^nd H2{x) in Fig. 
13 From Fig. ^ it is apparent that the high energy transfers are cut off because of the nuclear form 
factor at values lower than the limit imposed by the upper bound on the WIMP velocity, which 
increases with n. In the case of "^^{x), one clearly sees that the peak value decreases with n. Even 
though the area under the curve remains roughly independent of n, the portion available to direct 
detection decreases because of the nuclear form factor. In the case of H{x) one notices a change in 
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Figure 2: The function H[x) giving the effect of the velocity distribution on the modulated difi'erential rate. 
From left to right n = 1 and 2. Note the change in sign and the fact that the amplitude decreases with 
incensing n. 



sign. The low Q section tends to cancel the high Q part, when one integrates over Q to get the total 
event rate. Furthermore we see that the modulation amplitude is decreasing with n. By comparing 
H{x) and H2{x), one can see that the effect of the higher Fourier components in a is negligible. 



APPLICATIONS 



As we have already mentioned, the absolute rate depends critically on the specific nature of the 
WIMP, e.g. on the SUSY parameters in the case of the neutralino. It also depends on the structure 
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Figure 3: The function H2{x) giving a higher order effect on the modulated differential rate. From left to 
right n = 1 and 2. Note again the change in sign and the fact that the amplitude decreases with increasing 
n. 



of the nucleon. In the present work we will not be concerned with those very important aspects (see 
e.g. Refs HmSlilEi on how one deals with such issues). The event rate is proportional to the 
WIMP density in our vicinity, which is not modified by including the coupling between dark matter 
and dark energy as in our model. In any case, we will focus here on the aspects affected by the 
WIMP velocity distribution. 

The differential rate discussed in the previous section depends on the nucleus via its form factor and 
its mass. It also depends on the WIMP mass through the reduced mass Hr entering the parameter 
a. For our numerical study we will focus on ^^''l, which is one of the most popular targets employed. 
The nuclear form factor we use was obtained in the shell model description of the target and is 
shown in Fig. 0] 
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Figure 4: The form factor F'^{u) employed in our calculation, u is the energy transfer to the nucleus in units 
of Qo, i.e. u = Q IQo, with Qo = 64 keV. 



The part of the differential rate associated with indicated by dtcoh/du, is shown in Fig. I^lfor 
two WIMP masses = 30 and 100 GeV, and n = 1, 2. The explicit results of this figure confirm 
those derived by inspection of Fig. ^ 




Figure 5: The quantity dtcoh/du for n = 1 on the left and n = 2 on the right. From top to bottom = 30 
and = 100 GeV. 



The total (time averaged) rate is given by: 

tcoh = [ ^^^!^du, (37) 
J„^^^ du 

where Umin is determined by the detector threshold and Umax = {nyescY I ci^ by the maximmn WIMP 
velocity. By including both ^'o(a\/u) and H{a^Ju) wc can cast the rate in the form: 

Tcoh = tcoh (1 + hcoh cos a) 

hcoh = - — / du. (38) 

tcoh Ju„i„ du 

Integrating over the energy transfer, assuming either no detector cut off {umin = 0) or a cut off 
of Qth = 10 keV, we obtain the results shown in Fig. One can see from Fig. that, except 
for the case of light WIMPs, the total rate is decreasing with increasing n. The reason is that, 
as we have seen in the previous section, even though the total area under the curves of Fig. ^ is 
independent of n, the nuclear form factor dumps out the high u components. Similarly, the area 
under the curves of Fig. El decreases with n. As expected, for a given n the rate decreases as the 
energy cut off increases. From Fig. El we see that the modulation amplitude h rapidly decreases as 
n increases. This is expected in view of Fig. „ For a given n the modulation amplitude increases as 
the energy cut off increases. The reason is that h essentially is the ratio of the modulated amplitude 





divided by the unmodulated amplitude. Both decrease with increasing Qmin, but the denominator is 
decreasing much more rapidly. In other words the increase in the modulation occurs at the expense 
of the number of counts. 



CONCLUSIONS 



According to our present understanding of the evolution of the Universe, the main contribution 
to its energy content comes from two sectors, the dark energy and dark matter, that have not been 
directly observed. Within the majority of the models that have been proposed for the description 
of these sectors, there is no coupling between them other than the gravitational one. In general, 
this very strong assumption is not supported by some reasoning based on fundamental properties 
of the model, such as symmetries. On the other hand, a coupling between the two sectors may 
even be desirable, as it may provide an explanation of the coincidence problem, i.e. the comparable 
contributions of the two sectors to the energy density of the Universe today. 

The main motivation for this paper has been the wish to explore the direct observational con- 
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Figure 7: The quantity hcoh for Qmin = at the top and Qmin = 10 keV at the bottom. From left to right 
n = 1 and 2. 



sequences of such a coupling. We modeled the interaction between dark matter and dark energy 
by assuming that the mass of the dark matter particles depends on the scalar field whose potential 
provides the dark energy. The fact that the dark energy field, if it plays a dynamical role in the 
cosmological evolution today, must be effectively massless at length scales below the horizon means 
that its presence results in a long-range attractive force in the dark matter sector. This can have sig- 
nificant implications for the mechanisms of structure formation. The effect that has been of interest 
to us is the modification of the standard isothermal Maxwell-Boltzmann distribution of dark matter 
in the galaxy halos [s^ Is^ . The main modification is that the characteristic dark matter velocity 
can be increased significantly in the presence of the additional force. As the velocity affects directly 
the detection rates of the various experiments that search for dark matter, a detailed calculation of 
these rates, taking into account the new interaction, is important. 

The modification of the velocity distribution has consequences for direct WIMP detection. Re- 
garding the (time averaged) event rates, our results depend on the WIMP mass. For light WIMPs 
we find an increase of the rates by about 50%, independently of the detector energy cutoff, if the 
new force is stronger than the gravitational force by a factor n = \/3- For larger masses, however. 



the new force leads to a substantial decrease in the rates. The reason for this is that large energy 
transfers are inhibited by the nuclear form factor. The consequences of the new interaction are more 
pronounced in the case of the modulation amplitude. For light WIMP masses the modulation is 
decreased by an order of magnitude for an ideal detector (zero energy threshold) for k = VS relative 
to K = 0. For heavy WIMPs the decrease is about a factor of 4. The effect persists, but is somewhat 
less pronounced, in the case of a detector with a finite energy threshold, e.g. about 10 keV. 

We should emphasize that the average WIMP energy, in addition to its linear increase with the 
WIMP mass, also increases proportionally to 1 + (see Eq. Q). The same holds for the maximum 
WIMP energy, which depends on the escape velocity and scales with the same factor (see Eq. (fl^ 'l. 
This provides the opportunity of planning novel experiments other than those involving nuclear 
recoil, e.g. experiments detecting transitions to excited nuclear states in the MeV region. The 
relevant rate may be enhanced significantly due to the tail of the velocity distribution. 
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